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Abstract. Solutions of the stationary axisymmetric Einstein equations describing 
the interior of circularly rotating dust are investigated in order to study their potential 
applicability as galaxy models. It is shown that such interior solutions cannot be used 
as global metrics without becoming unphysical in certain regions of space. Although 
definite results concerning the non-existence of a continuation into a vacuum can only 
be derived for interior solutions describing rigidly rotating dust, the present analysis 
exhibits that the case of non-rigidly rotating dust would in general also be inadequate 
as a physically reasonable galaxy model. 
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1. Introduction 

There has been some recent interest in a general relativistic model of Cooperstock & 
Tieu fr, which attempts to fit rotation curves of spiral galaxies without invoking dark 
matter. This model is constructed such that it gives an excellent fit to observed rotation 
curves. However, there have been serious concerns about the consistency of the model 
in general Ej . 

By using some basic results presented in where solutions of the Einstein 
equations describing axisymmetric rotating dust were calculated, the applicability of 
the corresponding solutions to serve as galaxy models is studied in this paper. Such 
axisymmetric solutions have been used as toy models for galaxies during the last years, 
e.g. in pQ and [H] , in order to investigate whether general relativity can give new insights 
into the problem of the distribution of dark matter in our universe. 

The crucial question concerning the usefulness of axisymmetric dust solutions is 
whether the assumption of galaxies being stationary axisymmetric objects is in some 
way realistic. Since a galaxy is a dynamic object and a spiral structure is surely not 
invariant under arbitrary rotations, the above seems barely reasonable. Furthermore, 
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circular flow of dust makes only sense for the galactic disk, but far off the galactic plane 
such dynamics would become rather unphysical. This has already been pointed out in 
[Hj, where it was also shown, that the models in PI Ej indeed describe rigid rotation. 
Nevertheless, for a slowly rotating galactic disk where details of the spiral structure 
can be 'washed' out, i.e. where variations in the time coordinate t and the angular 
coordinate <p are rather small in comparison to variations in the other coordinates, the 
assumptions described above might yield useful approximations. A restriction to the 
galactic disk may also be reasonable in connection with the debate about dark matter, 
since the main question is then whether an inner galactic disk solution can be continued 
into a vacuum without the need for a dark matter halo. Indeed it will be show below 
that solutions describing simple axisymmetric configurations of rigidly rotating dust 
cannot be continued into a vacuum solution without involving singularities. 

Furthermore, also axisymmetric dust solutions with non-rigid rotation will be 
considered. In comparison to rigid rotation models, these solutions have a further degree 
of freedom expressible by an arbitrary function of one variable. A general answer about 
transitions into vacuum solutions will not be given, but according to what is currently 
known about vacuum solutions for Dirichlet boundary data [HI Ej, such transitions 
will in general exhibit singular behavior on the boundary. 

2. Einstein equations for static axisymmetric metrics 

Throughout this paper, geometric units are used such that 8nG — c — 1. In the case of 
a static axisymmetric problem as treated in the following, the metric can be assumed 
to have the simplified form 



where the functions W, V, R and N depend only on the coordinates /x and v. Details 
about this simplification can be found in, e.g. jU]. The form of the metric (0) includes 
the assumption that dt is timelike and that the planes orthogonal to dt and d v , i.e. 
and d p , are spacelike. This obviously excludes the possibility that the metric can have 
highly gravitating domains like ergoregions. 

In this paper, calculations are performed in the tetrad 



E = e~ w d t , E 1 = e w - y d,, E 2 = e w ~ v d u , E 3 = —[Nd t + d v ] , (2) 



where g(Ei, Ej) has the form of the Minkowski tensor. This semi-orthogonal basis is 
used since it allows some simplifications in the derivation of the curvature forms and 
the Einstein equations (see also 10J ) . However, the only apparent difference to the 
coordinate basis is a slightly different representation of the used tensors. 

The stress energy tensor 

The energy tensor of circularly rotating dust (i.e. a pressureless fluid) is defined 

via 




(1) 





(3) 
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where the timelike 4-velocity of the fluid Z is given by 
Z 



cosh — En — sinn — £3 
2 2 



(4) 



and p, as well as A, depend only on p and v. Here, A has been introduced as a degree 
of freedom for the 4-velocity Z, but if En and -R-E3 are interpreted as some kind of time 
and angular vector, respectively, then the value 



1 , A 
X = — tanh — 
R 2 



(5) 



can be interpreted as the angular velocity of the dust cloud in the frame . 

Furthermore, when writing Z as parallel to df — Qd^, where Q is the angular velocity 
seen by an observer who is not rotating with respect to the background of fixed stars, 
Q and X are related via 

(6) 



X e 2W (l + {lN) 

which can also be written as 

X 



n 



-2W 



XN 



(7) 



Now, the condition, that d t — ttd v must be timelike, which is equivalent to demand Z 
being timelike, can be written as 

n 2 R 2 < (1 + {lN) 2 e 4W . (8) 

According to ©, this condition is equivalent to \RX\ < 1. This inequality has the 
relatively simple physical interpretation, namely that the velocity XR of a particle at 
radius R in the frame (J2J) cannot exceed the speed of light. 

As a consequence of Bianchi's identity, the tensor T must obey the condition 
divT = 0. Since p only depends on the coordinates p and u, this is equivalent to 
the equations 



p div Z 
pVzZ 








(9 a) 
(96) 



where V denotes the Levi-Civita connection on the considered Lorentzian manifold. 
Due to the form of the metric, (}9"aj) is automatically fulfilled, whereas (|9^jl implies the 
further conditions 



pe 



w-v 



AT D 

cosh A + e 2W -£ sinh A - (cosh A - l)-£ 
2R 2R 



w-v 



W v cosh A + e 



^^sinhA-(coshA-l)^ 











(10a) 
(106) 



The indices denote partial derivatives. In the regions where p 7^ 0, it will be 
advantageous to use these equations instead of a subset of the Einstein equations. 
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Einstein equations 

With the stress energy tensor according to ©, the Einstein equations can be 
represented as 



a 2(w-y) -^W 







Xw-v) 



£{W-V)\ 



J.(W-V) 



a 2(W-V) 



a 2(w-v) 



+ 2W„ 



+ 2 



+ e' 



i\T 2 + TV 2 ' 



J2W 



R 



+ N w + 4W„N„ + AW V N V - 



R 2 

R^N„ + R y N v 



v m + v w + wi + w;, 



AW > 



AT 2 + N 2 



4i? 2 



p sinh 2 



p cosh A 
= —p sinh A 



A 



W 2 -W 2 - ^ 
" " R 



U V U A " I 'I' 



R 



R 



2R 2R 



N 2 - N 2 

V " 



i? 



i? 



2i? 2 



4i? 2 
0. 



[lla) 
(116) 

(11c) 

[lid) 

(lie) 

(11/) 



If i? 2 + R 2 7^ 0, which makes sense if R is interpreted as some kind of radial distance 
from the rotation axis, then (jllal e. /) can be used to find explicit expressions for and 
V u as follows 

1 _ _ R U W 2 - R U W 2 - 2R V W„W U 

v,=-dMRl + Rl) + — ' " 



4R(R 2 + Rl) 



+ e 



m R,N 2 u -R,Nl-2R v N fl N u 



12a) 



4R(R 2 + Rl) 



V u =-d u \n(Rl + Rl) + 



+ e 



AW 



, R U W 2 + R U W 2 - 2R P W P W V 

i + Ry) 

—R U N 2 + R v Nl - 2R tl N tl N u 



4R(R 2 + Rl) 



126) 



AR{Rl + Rl) 

The integrability condition for §12 al b) is a consequence of (jlUal fe) and (JTT^-c). 
Furthermore, (jlldj) can be derived from those equations. Thus, the original problem of 
solving (jlld - / ) can be reduced to finding solutions of f|10a| 6) and (jllal -c). 
By introducing the flat pseudo-metric 

gP = dfi 2 + dp 2 (13) 

and defining the derivative operators div^, and A v as well as the scalar product 
(•, •) like in the common case of two-dimensional cartesian coordinates, these equations 



can also be rewritten as 
P 

A r R = 



, . ,-vp TT r 9W sinh A T cosh A — 
cosh AV r W + e 2W — V V N — V v R 



2R 



2R 



(14a) 
(146) 
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-idiv p [2RV V W] + — (V V N,V V N) = e^-^pcoshA 



div 



v 



J2W 



R 



-V V N 



R 2 
2e 2W 



+ ^— (V V W, V V N) = -e 2 ( y ~% sinh A . 



(14c) 
(lAd) 



It should be emphasized, that similar equations could be derived for a fluid with pressure. 
Furthermore, these equations would also be valid in the vacuum case - however, (jl4aj) 
would then be trivially fulfilled. 

In the upcoming calculations of the interior solution, the fact that (jl4aj) is not 
trivial plays an important role. 



3. Solutions of the Einstein equations 

When writing p = A V U for some Potential U, the structure of equations (jMol oO is 
invariant with respect to coordinate transformations = z/) where £ and rj 

are harmonic functions satisfying £^ = r] u and £ v = — r^ M . In other words, coordinate 
transformations that can be written as ^ + ir] = f(fi + iv) for some holomorphic function 
/■ 

A helpful consequence of this fact is that R, together with a function Q obeying 
Qfj, = R v and Q u = —R^, can also be considered as coordinates, since by (jl4fej) R is 
harmonic in p and v. Vice versa, R can be taken as any harmonic function in p and 
v, since a choice of R is nothing else but a choice of the parametrization. The only 
restriction is that R should remain positive, when this function is considered as a radial 
distance from the rotation clXIS, clS it will be done in the following. This interpretation 
of R is quite natural, because in the case W = 0, N = the metric looks like the flat 
Minkowski metric in the coordinates (x, y,z) = (R((jl, v) cosy?, R(/jl, v) sin if, Q(p, v)). In 
fact, this is not a real restriction to the choice of R, but merely a restriction on the range 
of the values of the coordinates p and v. By this interpretation of R, the associated 
function Q can be interpreted as the distance from the hypersurface defined by Q = 0, 
that can be set as the symmetry plane - if the solution is assumed to have a reflection 
symmetry. This makes R and Q interpretable as some kind of cylindrical coordinates, 
called Weyl's canonical coordinates. 



3.1. Rigid rotation and static solutions 

For Q = 0, which due to (0) and (0) is equivalent to A = 0, the dust is not rotating 
with respect to the background of fixed stars and can thus be called static. This section 
is dedicated to the relation between solutions describing rigid rotation and stationary 
solutions by the use of co-moving coordinates. This transformation makes (|14al -c0 look 
less complicated for A = 0. This procedure to find solutions has also been used in 
PU 121 E] . Furthermore, it will be shown that solutions describing rigid rotation obtained 
from static solutions by this procedure will lead to unphysical metrics unless the metric 
functions obey certain inequalities. 
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3.1.1. About co-moving coordinates Solutions of (|14oI -gO describing rigid rotation are 
obtained by setting the 4- velocity Z parallel to d t — Sld^ where fl is constant, but 
nonzero. Since matter is assumed to move on timelike geodesies, condition (jHJ) must be 
satisfied where p ^ 0. However, for any asymptotically flat solution the right hand side 
of (JBI) should go towards 1 for large values of R or Q, whereas the left hand side can 
be arbitrarily large. This implies that the rotating dust cloud must be contained in a 
domain with R < R max . Note that the conclusions above also hold when pressure is 
included. 

The co- moving coordinates are chosen in a way to have c\ = dt — fld^. The simplest 
way to accomplish this is to set (p := p — fit and to take the identity mapping for the 
other coordinates. Then also Z = E and the transformed metric looks like 

g = [ e ^(i + fiiv) 2 - e - 2W Q 2 R 2 ] d? - 2 [e 2W (l + QN)N - e~ 2W QR 2 ] d* dp 

_ {e -W R 2 _ e 2W N 2 )d ^2 _ e 2(V-W) {d ~2 _ ^-2) ( 15 ) 

Thus, where (JHJ) holds, c\ is timelike and the transformed metric (fl3j) can be written in 
the same form as with the metric functions 

e 2W = e 2W {1 + QN) 2 _ Q -2W Q 2 R 2 (1 g fl) 

~ _ e 2W (l + flN)N -e~ 2W flR 2 

~ e 2W (l + QN) 2 - e~ 2W fl 2 R 2 ( ' 

e 2V 

e 2y = - (16c) 

e 2W (l + QN) 2 -e~ 2W Q 2 R 2 y ' 

R 2 =R 2 . (lQd) 

Hence, by the use of co-moving coordinates all solutions describing rigidly rotating dust 
can be reduced to the stationary case A = 0. The transformation back into non co- 
moving coordinates can be obtained from the above by simply replacing fl with —fl 

as 

Q 2W = Q 2W {1 _ Q ^2 _ Q -2W Q 2 R 2 

N e 2 ^(l-flN)N + e- 2 ^flR 2 
e 2W {l - flN) 2 - e~ 2W fl 2 R 2 

e 2V 

e 2V = J (17c) 

e 2W {l-QN) 2 -e~ 2W Q 2 R 2 

R 2 =R 2 . (17 d) 

It can be read directly from (jlfiqjl that e 2W > is equivalent to (jHJ) and similarly, 
according to ()17a|) . e 2W > is equivalent to (JHJ) for the metric functions l\liia\ -d) with 
— fl instead of fl. Hence, a stationary solution satisfying 

< e 2W (18a) 
< e 2 ^(l - flN) 2 - fl 2 R 2 e~™ (18b) 

can be transformed into an interior solution describing rigidly rotating matter with 
timelike d t and Z. Of course, by what was said at the beginning of this paragraph, such 
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a solution would not be asymptotically flat or make no sense as a global solution from 
a physical point of view. It was already shown in [2 j that globally defined rigid rotation 
solutions are plagued by unphysical properties. Further, it might also be convenient to 
demand that SIN and SIN should be of positive sign in order to prevent the solution 
from containing some counter-rotating sources of matter or singularities. This leads to 
the additional conditions 



< SIN 

„2W 



1 _ UN) 2 - e~ 2W Sl 2 R 2 < e 2W (l - SIN) . 



(19a) 
(196) 



A further condition for a reasonable solution would also be to demand that should 
be spacelike, which is equivalent to demand d$ being spacelike, i.e. 



e~ 2W R 2 



e 2W N 2 >0. 



(20) 



3.1.2. Einstein's equations for static dust In the static case, and for p ^ 0, (|14ajl 
requires setting W equal to a constant. Since the addition of an arbitrary constant to 
W would just result in a scale transformation, this constant can be set equal to zero 
without any loss of generality. Thus, it can be assumed W = henceforth. With these 
simplifications ^V2a\ b) can be written as 



-d, ln(^ + R u ) + - — 



l -d v \n(Rl + Rl) 



-RuN 2 



R„N 2 



2R,N,N U 



AR{R? + Rl) 



Furthermore, (|14d c?) reduce to 
V V N 2 



R 
div r 



2V 

e p 



V F N 
R 



0. 



(21a) 
(216) 

(22 a) 
(226) 



Obviously, these equations are completely determined by the function N, which itself 
must obey (|226j) . An important property of these equations is that, according to (|22aj) . 
the energy density p is always positive. In ()14c|) this cannot be assumed a priori - and 
as calculations will show, in the case of non-rigid rotation the mass density can indeed 
become negative. 

Some examples, how explicit solutions of that equation can found in some specific 
parameterizations are given in appendix |Appendix B 
solutions of (12261) is mentioned in 13.1.31 



Another way to construct 



3.1.3. Rigidly rotating dust as a galaxy model This section about rigidly rotating dust 
solutions will be concluded by a short discussion about the applicability of those metrics 
as galaxy models, as it has been done in [HE]. Since in both papers an attempt is made 
to give global solutions, they surely must, as mentioned before, contain domains where 
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the metric becomes unphysical. Where the problem lies in |5| will briefly be analyzed 
by the use of KTM b) and fiTM b). In [TJ |3] it has been set W = 0. Hence, (IT£al) holds, 
so according to (fT8ljl and flffl fr), fiiV must be bounded by < O/V < 1. Thus, (fTslj) 
will not hold if R > R max for some R ma x- This will happen for any choice of Q ^ 0. As 
a result, the metric in non co- moving coordinates has a time-coordinate which becomes 
spacelike as soon as R > R ma x- 

To avoid this problem, the metric for rigidly rotating dust must be restricted to a 
domain with R < R max where fll8a( fe). (jl9a| fe) and hold, and continued by another 
metric. Remembering that (jl9aj) is only necessary for rigidly rotating solutions, the 
continued metric could also be a metric describing rotating dust, but with non-rigid 
rotation and an angular velocity that falls off to zero for large values of R and Q. 
However, any transition into such an axisymmetric dust solution would still contain 
another source of unphysical behavior, namely that circular orbits are still geodesies 
according to (pJfe)) . This is only physically reasonable in the vicinity of the galactic 
plane Based on this observation, a further restriction of the solution to a domain 
with \Q\ < Q max , when Q is interpreted as the distance from the galactic plane, and 
a transition into a vacuum state appears to be the simplest choice to obtain physically 
reasonable solutions. Unfortunately, as will be shown later, this attempt will result in 
some other singularities. 

It has already been shown by Van Stockum ^T] that interior solutions for rigidly 
rotating dust can be obtained by axisymmetric solutions of Laplace's equation in three 
dimensions. This relation does not only allow more general methods to find solutions of 
(122 6|) than shown in |Appendix B[ but can also give some deeper insight into the behavior 
of the solution due to corresponding theorems about harmonic functions. The relation 
mentioned above can be constructed as follows. Using the linear mapping J : M 2 — > M 2 
describing a rotation of |, i.e. J(x,y) = (—y,x), an immediate consequence of (j22fc|) is 
that, according to Poincare's lemma, a function H can locally be defined by 



up to an arbitrary constant. Especially when N is defined for /z, v G U where U is a 
simply connected domain, H can also be defined on whole U. This function H must 
then satisfy 



which is nothing else but Laplace's equation for an axisymmetric function in three 
dimensions in the coordinates (x,y,z) = u) cos (p, R(fJ,, v) sin <p, Q(fi, u)). Two 

quite fundamental consequences of the relation (|2^j) are listed below. 

One of the most important observations is a direct consequence of the fact that 
harmonic functions have no isolated extrema. Thus, any non-trivial globally defined 
asymptotically flat solution iV of (|22frj) must contain a singularity - if AT is continuous, 
then its first or second derivatives must contain some incontinuity. This is not really 
astonishing, since such a solution would in fact describe a non-rotating dust cloud 



V V H = -JV V N 



(23) 



div v [RV V H] =0, 



(24) 
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without any inner pressure to keep the particles in their place. Hence, from a physical 
point of view, the logical consequence is that there must be some kind of singularity with 
a negative mass, as shown in |2j, to prevent the cloud from collapsing. Furthermore, 
the presented discussion illustrates that it would be a hopeless undertaking to model a 
spheroidal galaxy of type EO in the Hubble sequence as gravitationally stabilized dust. 

Another consequence is that the whole spacetime cannot be covered singularity- 
free by restricting the rotating dust solution to an axisymmetric region where \x and 
v live on a domain on which Poincare's lemma can be applied. Due to ()22aj) . such a 
transition would yield N = const, and thus H = const, on the border of the domain 
containing the dust cloud. As a consequence, H, and therefore N, must be constant on 
the whole domain, implying that there is also p = 0, i.e. a vacuum, inside the domain 
according to (j22aj) . For the use as galaxy models, it would probably be sufficient to just 
consider interior solutions on simply connected domains. Nevertheless, the nonexistence 
of transitions into a vacuum for more general measurable domains can also be verified 
under some further assumptions as follows. Let /i, v e U, where U is bounded with a 
piecewise smooth Lipschitz boundary dU = Ay + T, where Ay is the intersection of U 
with the rotation axis and T the 'real' boundary of the domain containing the rigidly 
rotating dust cloud. Furthermore, N/R and V V N/R should be well-defined as Sobolev 
H 1 functions on U and (V V N, n) = on Ay, where n is the unit normal direction. These 
assumptions are not too devious, since they guarantee that the metric does not become 
singular on the rotation axis (if W satisfies similar conditions) and also, if the ratio 
N/R is small enough, that d v remains spacelike. By using (|226J) and Gauss' theorem, 
it follows 



A direct consequence of demanding V V N = on T is V^A" = and therefore a vacuum 
inside U. 

As the calculations above show, a cloud of rigidly rotating dust on a bounded 
domain cannot exist without any forces stabilizing that matter configuration. 
Furthermore, some embedded vacuum solutions, even including exotic objects like 
toroidally shaped black holes, would not do the trick, since any continuation into a 
vacuum solution would be insufficient. Thus, only continuations into other nonempty 
solutions, like e.g. a dark matter halo, are possible choices. Hence, claiming that 




(25) 
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solutions considered in this section are realistic models for the galactic disk could rather 
be viewed as an argument in favor of the existence of dark matter than against it. 

3.2. Solutions describing arbitrary rotation 

In the following it is presupposed that A 7^ in all points. The corresponding interior 
solutions for non-rigidly rotating dust have already been elaborated by Winicour in 
1975 jl], a derivation adapted to the form of the metric used in this paper is given 
Appendix A The solutions are completely characterized by an arbitrary positive 



111 



function of one variable q and a solution of (J22fcj) Furthermore, W and A" can be 
expressed as functions of R and X © as 

e~ 2W = T-^SW (26) 



1 - R 2 X 2 

With constants Xo and Nq. The value of X is related to q and ^> by 

M> = R 2 f(X) + h(X) , (28) 

where 

f\X) = 2-X*h'(X), h'(X) = ^- y (29) 

Since 1 and R 2 are solutions of ([22 6[) . it is not a real obstacle that, according to (f2~9*]). 
/ and h are only defined up to some constants. This could e.g. be fixed by demanding 
f(X ) = and h(X ) = 0. The important question is now, whether (|25|) can be solved 
for X. Assuming 

R 2 f'(X ) + h'(X ) = 2R 2 Q + (1 - R 2 X 2 )h'(X ) + (30) 

for some value Rq, the theorem of implicit functions guarantees that (|2*8*j) can really be 
solved for X as long as \I> and \R — Rq\ are small enough, since for X = Xq the solution 
of (J2HJ) is given by ^ = . Based on the theorem, the gradient V v X can be expressed 
as 

V P X = * [V"* - 2f{X)RV v R] , (31) 



with x as defined in (|A.9j) . Because the solution ^ = for X = X is given for any 
choice for R, (|28j) can be in general only be solved for all values of R if ^ satisfies some 
condition C(R) > \^>\ for a positive function C depending on q. Nevertheless, in the 
special case where h' is positive with h'(X ) > and X 2 h'(X) < 5 < 2, condition (j371jl 
would always be true and (|2^j) can be solved globally without any further restrictions 
on \I>. 

Furthermore, the energy density can be expressed as 
e 2(v-w )p = Q'(X) 2 - R* [2XqjX) - X 2 q\X)f 



Aq 2 R 2 X 2 
AR 2 l V X l 



(32) 
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By means of ()31j) . p can also be written as 



e Kv-w )p 



(1 + R 2 X 2 )q\X) - 2R 2 Xq{X) 
(1 - R 2 X 2 )q'(X) + 2R 2 Xq(X) 

(1 + R 2 X 2 )h'(X) - 2R 2 

(1 - R 2 X 2 )h'{X) + 2R 2 2R 



2R 



- f(X)V v R 



- f(X)V r R 



(33) 



Since X and R only depend on the variables \i and i/, p can also easily be expressed by 
using the gradient assigned to the metric g instead of V v as 

h'(xy - ^VW^ adX)gradX) 

(34) 



4i? 2 



1 



AR 2 



g (grad ) + # 2 grad f(X), grad fi(Jf) - i? 2 grad /(X)) 



4. Discussion 



The function q has been introduced as a degree of freedom. Since q depends on the 
angular momentum X, and W is constant in the stationary case X = 0, the function q 
must contain some information about the rotation law of the dust cloud. This can be 
seen in a more formal way as follows. According to (J7|) as well as (126)) and (|27|). the 
angular velocity Q obeys the relation 

1 _ q(X) 1 
Q X 2 



X A *M V 
dy N 

Xq 



X 2J y y 
Taking the gradient of this equation yields 
1 



(35) 



grad 



n 



X 2 2X 



gradX = 1 [2g(X) - Xg'(X)] grad 1 



(36) 



Thus, the value 2g(X) —Xq , (X) can be taken as an indicator how strongly the rotation 
of the matter deviates from rigid rotation. Especially if that factor vanishes the rotation 
is rigid - of course, there would also be no non-rigidly rotating matter if gradX = 0, 
but because of there would be no matter at all. 

As aforementioned in 13.1) the energy density can be negative if W is not constant. 
According to the energy density is positive only as long as 

h\X) 2 - R 4 f'(X) 2 > . (37) 

This relation can also be written as 



R 2 < 



h'(X) 



2 - X 2 h'(X) 



(38) 



X[2q-Xq'(X)] 

provided the denominator is different from zero. The cases where ()38)1 is not well 
defined are X = and 2q — Xq'(X) = 0, i.e. when there is no or just rigid rotation - 
presupposing q' as continuous. Hence, if the rotation is differential, the deviation from 
rigid rotation implies a restriction on the maximal radius of the rotating dust cloud, 
which is inversely proportional to the deviation. This inequality above is a further 
restriction for R besides \XR\ < 1, since for arbitrary h' the latter cannot be derived 
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from (jHZJ) or ()38|). Especially, when > R 2 > \h'(X)\, p will be negative. Hence, all 
physically reasonable solutions must obey R 2 < \h'(X)\. 

However, when assuming < X 2 h'(X) < 1 then ()38j) also implies \RX\ < 1. This 
implication is strict in a sense that for any X with X 2 /i'(X) > 1 there exists an R 

RX 



that obeys (|5Sj) but i?X > 1. A direct consequence of the discussion above is that all 
matter must be contained in a domain with a finite radius when h! is bounded, as in 
the case of rigid rotation. 

The restriction of the dust cloud in the radial R direction will also follow if the 
solution converges smoothly enough to the asymptotic flat space. Since for R — > oo it 
is necessary to have X/Q — > 1, implying by that q — > with g^X) = 1 — (No/ 
X )X 2 for X 7^ and q^ = 1 otherwise, because X = implies N = when 
asymptotically X — > 0. According to (jHHjl R 2 must then be bounded by |X /X | or 
respectively, which is a contradiction to R — > oo. Thus, the matter cannot be contained 
in a region with arbitrarily large values for R. 

It was mentioned that 

(1 + R 2 X 2 )h'(X) - 2R 2 = h'(X) - [2 - X 2 h'(X)]R 2 = (39) 

can have nontrivial solutions. These mark the points of the spacetime where p is 
vanishing, assuming j^\7 v ty is not singular in these points. Furthermore, according 
to (jUID), the equation 

(1 - R 2 X 2 )h'(X) + 2R 2 = h'(X) + [2 - X 2 ti(X)}R 2 = (40) 

marks the points where the energy density is potentially singular. Regarding ()30|) . then 
(pfUj) marks the values (R, X) where the implicit function theorem cannot be applied. 
This is because V V X would become singular and additionally the implicit equation 
might possibly have a branch-point there. 

If R > 0, the energy density cannot be positive for h'(X) = for obvious reasons. 
Thus it can be assumed that h' will not change the sign on the whole domain with 
p > 0. To avoid singularities and being able to find coordinates, where p can vanish to 
zero, it is (according to (J3T?|) and (jUIlO quite advantageous to choose h! as positive on 
that domain. Remembering that R and X are both functions of the two variables p and 
u, and that (|14aj) is only non-trivial for p ^ 0, (J39)) can be seen as the specification for 
the border of the region where the interior solution makes physically sense. By choosing 
q, or h! respectively, and ^ in a suitable way, interior metrics can be found, where p 
indeed vanishes on the surface of a bounded domain containing the dust, which could be 
seen as the edge of the galactic disk. This situation allows a somewhat more formalistic 
approach to the dark matter problem in terms of general relativity, namely whether a 
smooth transition from the rotating dust into an asymptotically flat vacuum solution can 
be found or does not exist. Though some axisymmetric vacuum solutions are known [12J, 
especially including solutions outside an infinitesimal disk EH Hi] , which could yield 
good approximations for the vacuum solution outside a galactic disc (should it exist) 
there does not yet seem to be a theorem that could guarantee either the possibility or 
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the impossibility of a smooth transition. However, what has been exhibited by recent 
developments [HJ IE] is that the asymptotically flat vacuum solution outside a bounded 
region for given Dirichlet boundary data is unique and the existence of a solution only 
seems to be guaranteed, if these data satisfy certain conditions. Furthermore, in [7] an 
approximation scheme was presented, how solutions can be obtained from it's Dirichlet 
boundary data. Since smooth transitions not only require a correspondence of the 
Dirichlet data at the boundary, but also of the Neumann boundary data, the problem of 
finding an attached vacuum solution for an arbitrary choice of q and ^ is overdetermined 
and would thus in general create some kind of mass shell at the boundary, as it also was 
mentioned in |Hj. 

4-1. The limit of slow rotation 

In the case of slow rotation where X is small, the solution in 13.21 should in some way be 
well approximated by a static solution, as described in 13.1.21 Thus, a solution with an 
arbitrary rotation law will be expanded in X to lowest order. 

In the following it is assumed that X = and q(0) — 1. In fact, q(0) can be taken 
as any nonzero number, since multiplication of q with some nonzero number would 
simply result in a scale transformation. The expansion of (|26J) and (|2"Yj) then yields 

e 2W = 1 + 0(X 2 ) (41) 



and 



iV = | + iV + O(X 2 ), (42) 



where f)28|) has been used. 

Furthermore, an expansion of (|33|) will unveil the result 



9V h'(Q) - 2R 

e 2V p - 



h'(0) + 2R 2 



V 7 ^ 



2R 



2 



O(X) . (43) 



Thus, when the rotation is rather slow and R is small in comparison to h'(0) the metric 
will indeed decay to a static solution as obtained in 13.1.21 Even further, (|43|) also shows 

that \ ^p- can be taken an indicator for the radius of the dust cloud. 



4-2. Back to rigid rotation 

Though the method of using co-moving coordinates described in 13.1.11 might be the 
easier attempt to find solutions describing rigid rotation, those solutions must also be 
included in those from 13.21 - in non co- moving coordinates. According to IJflfij) . rotation 
is rigid for q(X) = ^r, where an arbitrary factor has again been omitted. 

By fixing the factors X and iV according to (j3*3|) . it is now straightforward to 
calculate 

e 2 W= &__ n2R2 (44) 

N = — . (45) 
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Furthermore, (|28|) takes the form 

m = 

and p is given by 



4 



e *(v-w) p 



2R 



(46) 



(47) 



When comparing these results to the transformation rules (jl7al -<f) with W = and 
Y = ^(1 — fiiV), the metric is indeed the solution for rigid rotation described in 13.11 



5. Conclusions 



Besides the restriction that physically realistic solutions describing axisymmetric 
rotating dust should be restricted in the Q direction because, as already stated in 
the assumption of circular flow is unrealistic far off the galactic plane, the calculations 
have shown that it is also required to restrict the solution in the R direction. This 
unmasks the solutions in QjJ |3] once again as unphysical, though they seem to work 
quite well for the interior of the galactic disk. However, according to the expansions in 
14.11 the latter is not accidental, but simply because those metrics approximate non-rigid 
rotation solutions to first order. 

Furthermore, solutions describing rigid rotation cannot simply be matched with 
exterior vacuum solutions without becoming singular, making them a quite useless mean 
to claim the non-existence of dark matter, besides the fact that rigid rotation might 
possibly be quite inappropriate to describe the galactic rotation law, even though the 
measured rotation curves look like differential rotation due to the relativistic dragging 
effect observed by rotating objects. 

Even when non-rigid rotation of matter is allowed as in 13.21 things do not really 
improve in an optimal way. Though, in comparison to the case of rigid rotation, the 
interior solution can now indeed have an energy density p that vanishes on a compact 
surface, the resulting problem to find a matching vacuum solution is over-determined. 
Thus, the resulting metric would in general contain an infinitesimal mass shell enveloping 
the galactic disk. Nevertheless, it cannot be excluded at the current time, that smooth 
transitions could exist. Of course, also transitions into other solutions, especially such 
with more degrees of freedom should be analyzed in the near future, e.g. a dust solution 
where the 4- velocity Z is not restricted to be a linear combination of just dt and d^. 
But still, already this assumption would require that metrics of a more general form 
than have to be analyzed, making the chance to obtain some reasonable results that 
would finally conclude the debate about the existence of dark matter, quite difficult if 
not impossible at the current moment. 

Nevertheless, the results exhibit that galaxies should be treated as much more 
complex objects than just a circularly rotating cloud of dust and a final solution to 
match different solutions of the Einstein equations would also give some definite answers 
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concerning the existence of dark matter. Though this answer is still outstanding, the 
apparent incapability of solutions containing circularly rotating dust as the only source 
of gravity to serve as global metrics might give a further indication that more physics 
is needed to model galaxies, like pressure related phenomena, dark matter halos, effects 
caused by a massive central black hole or even non-axisymmetric dynamics. 



Appendix A. 



For p ^ 0, i.e. if matter is present, 
V V N and V V R as 



, , „i3 Trr ow sinh A -p „ _ cosh A 
cosh AV V W + e 2W — V V N 



implies a non-trivial relation between V V W, 
1. 



V V R = . 



Using the linear mapping J : 
is 

diy'Pj U AV7^T" owSinhA 



. . (A.l) 
2R 2R K ' 

R 2 introduced in 13.1.31 a direct consequence of (jA.l|) 



cosh A V W + e 



vPN _ coshA-1 " 



0. 



(A.2) 



2i? 2i? 

This equation guarantees that (|A.1|) can be solved for one of the values W, N or R, i.e. 
the integrability condition for V V W, V V N or V V R holds. 

Calculating the application of div v in ()A.2j) and using (jA.l)) to get rid of the term 
V V N leads to 



V V R V P A 



v, 



J 



V V R coshAV^A 



R sinhA'^ [' " 2R 
This equation is solved, when W is written as 

1 + cosh A 



2 sinh A 



0. 



W 



-bn 

2 



(A.3) 



(A.4) 



where q : 1R — > M>o is an arbitrary differentiable function and X is the function defined 
in (jSJ). As it will turn out, it is convenient to use the function X instead of A or Q. 
Hence, the following relations are useful 

2RX 



sinh A 



cosh A 



1 - R 2 X 2 
1 + R 2 X 2 



1 - R 2 X 2 ' 

Now (jA.l|) . (jA.4|) and (jA.51 6) can be used to calculate 



V P W = — 



2R 



R 2 X< 



1 - R 2 X 2 
2RX 



V V R + X (R,X) 



1 - R 2 X 2 
1 + R 2 X 2 



v v x 



l_ R 2 X 2 ^ ' ^1-^2^2 

where X) is defined by 

1 - R 2 X 2 )q'(X) + 2R 2 Xq(X) 



x(R,x) 

From this also follow fl2HJ) and (1271) . 



2Xq{X) 



(A.5) 
(A.6) 

(A.7) 
(A.8) 

(A.9) 
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The calculations above show that for any differentiable function q and a 
parametrization via R, the solution is determined by the function X, as the solution for 
rigid rotations was completely determined by a solution of (J226j) . Due to ()14c| g0. the 
function X cannot be chosen arbitrarily. The crucial point is now to derive an equation 
for X from ()14c| g0. As calculations will show, X is related to a solution of (J22fe)l . 

The procedure to derive the mentioned equation for X is as follows. First, from 
flMct cO a new equation can be obtained which does not contain the term p due to the 
the combination (jl4c)l sinhA + (|14cjj) cosh A = 0. This reads as 

= sinh A I - div v [2RV r W] + — < V V N, V V N) J 



+ cosh A ( div v 
sinh A 



J2W 



R 



-V V N 



2e 2W 

+ ^— (V V W, \7 V N) 
R 



R 



div v [2RV V W] +^^div v [e 2W V r N] 

R 



(A.10) 



R 2 
2e 2W 



+ sinh A^- (V V N, V V N) - cosh A 1 



R 2 



J2W 



V V R, 



R 



-V P N 



+ cosh A—— (V V W, V P N) . 
R 



When using (|A.1|) in order to reduce the number of terms, this can be rewritten as 




shlhA div* [2RV V W] + ^Adiv p [e 2W V v N] 

R 



R 



J2W 



(A.ll) 



Now, the terms cosh A and sinh A can be replaced by (jA.5|) and (jA.6|) yielding 







2X 



R(l - R 2 X 2 , 



div- [2RSTW] + „\ +R mL ^ [* 2W V V N] 



R(l - R 2 X 2 ) 



(A.12) 



By inserting fl2E|), $E3$) and (jA.Sjl . (jA.12j) can be written as 



2X 



R(l - R 2 X 2 



div 



v 



W V R- 



1 + R 2 X 2 
1 - R 2 X 2 



V r R-x(R,X) 



2X 



1 - R 2 X 2 



V V X 



1 + R 2 X 2 , 
+ -^tz — div 



v 



R(l 
1 

'R 2 
This equivalent to 



- R 2 X 2 ) 
V V R^ 



2RX 



R 2 X 2 



V y R + X (R,X) 



1 + R 2 X 2 
1 - R 2 X 2 



v v x 



(A.13) 



WX ^R + X (R,X)±±££;V' P X 



R 2 X 2 



R 2 X 2 



= -^[ X (R,X)V-X] + -(V-R, T - WT2 W 

~l + R 2 X 2 



1 + R 2 X 2 



v 



2RX 



R 2 X^ 



1 v 2RX 

rV R 'T^Wx' 



V 



1 - R 2 X 2 



1 

R 2 



V V R, 



2RX 



W V R) (A. 14) 



1 - R 2 X 2 



1 ^- 2 (V v R,x(R,X) i 



l + R 2 X 2 , 



R 2 X 2 



v v x 
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By straightforward calculation, this can be reduced to 



= div 



v 



R 



v v x 



+ 



2-X 



y v R,v p x) 



Since R is harmonic ()14fe|) . the equation above is equal to 
= div v 



X ^ X) V v X + f(X)V v R 



R 



(A.15) 



(A.16) 



where / is defined by 
be written as 



With the substitution (1281). what by use of (jA.9|) can also 



ds X (R,s), 



(A.17) 



when assuming f{X ) = and h(X ) = 0, ()A.16|) indeed turns out to be equivalent to 



div 



v 



R 



0. 



(A.18) 



This is the same equation as (J22foj) . 

By means of pijl. (|A.7|) and (jA.8|) can also be expressed as 



V V N -- 



X 



1 - R 2 X 2 



[f(X)-X}RV v R 



[l - R 2 X< 



( 



1 + R 2 X 2 
1 - R 2 X 2 



2RX - 



Rf(X) V H R 



1 + R 2 X 2 



+ 



(A.19) 



(A.20) 



1 - R 2 X 2 2 

The energy density can be obtained by inserting (jA.10|) into ()14c|) 



Appendix B. 

The differential equation (|225|) plays a central role when dealing with circularly rotating 
dust. Even though solutions can be obtained from axisymmetric harmonic functions 
in three dimensions, it would be quite practical to have some additional explicit 
expressions, how solutions can be found by separation of variables, when assuming that 
R is separable in \x and v. The so obtained solutions correspond to the most common 
coordinate systems that are suitable for axisymmetric functions. 

In the following some examples are given for R, Q and (j22fc|) in different coordinate 
systems. The designation of the systems refer to the parametrization in the case of flat 
space, i.e. when N = and W = 0. 



Cylindrical coordinates : 

R=u, Q = fM (B.l) 

i/* w + i/* w -* 1/ = (B.2) 
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Separating this equation with constant k 2 , k G C leads to 

u) = [A K e K » + B K e~^] [CWiM + D k vN x (kv)\ , (B.3) 

where J\ and iVi are the Bessel function of the first and second kind with order 1. When 
further assuming reflection symmetry, (|B.3|) can be simplified to 

v) = cosh(«/i) [C k vJi(kv) + D k vNx{ku)] . (B.4) 

Solutions of this kind have been used for the modelling of galactic rotation curves. In 
P it was claimed that such solutions are globally valid, producing different sources of 
unphysical behavior, mainly because rigidly rotating dust is described which must result 
in some unphysical behavior, as mentioned in j2] and discussed in 13.1.31 Nevertheless, 
as shown in 14.11 the results of Cooperstock & Tieu remain applicable in the vicinity of 
the galactic disc. 

Another way to produce solutions of ()B.2|) in cylindrical coordinates is in the form 
of the integral representation 

*(**,!/)= / dzC(z)- £ -j (B.5) 

Jr ((/j, — zy + v l ) 3 

with an arbitrary function of a complex variable C and a curve T in the complex plane. 

E.g. in |T| = M and an even function C was used to cover the whole spacetime. 

P araboloidal coordinates : 

R = fiu, Q = l(^-u 2 ) (B.6) 

/ii/* Wi + /iZ/tf M/ -z/tf M -/i¥„ = (B.7) 
Separation with the constant G C leads to the solution 

v) = fib> [A k Ji(k/j,) + B k Ni(k/j,)] [C k Ji(iku) + D k Nx(iku)] . (B-8) 
Reflection symmetry can be obtained by symmetrization like 

* K (n,is) =* K (n,i/) + * K (i/,n). (B.9) 

Spheroidal coordinates : 

R = R e^ sin v , Q = Roe 11 cos v (B. 10) 

sin v + sin v N vv — sin vN^ — cos vN v = (B-ll) 

Solutions of this equation with separation constant A(A — 1) G C, 3?A > —1 can be 
written as 

( [Aoe" + B ] [C cos u + D ] A = , 

* A (^, v ) = I sin 2 v [A x e {x+1) » + ^e"^ * (B.12) 

{ [C A P{(cos i/) + D x Q' x (cos v)] A ^ , 

where P\ and Q\ are Legendre functions of first and second kind. For arbitrary A, 
P\ and Q\ contain a logarithmic singularity at —1 and 1 respectively. Exceptions are 
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A = I G N, where p is a polynomial of degree I. When further assuming reflection 
symmetry, (jB.12|) reduces to 

A e^ + B A = , 

*x<jl,v) = { sin 2 ^ +1 (cos I /)[A 2m e( 2i+2 ^ (B.13) 
+B 2l+1 e^ 2l+1 ^] A = 2/ + l . 

Prolate spheroidal coordinates : 

R = Ro sinh /x sin v , Q = R cosh /x cos ^ (B. 14) 

sinh /x sin z/ + sinh /i sin v N vv — cosh /x sin vN^ — sinh /x cos vN v = (B.15) 

Separating again with A(A — 1) 6 C, 3?A > — 1 leads to the solution 

{[A cosh /x + .Bo] [C cos z/ + D ] A = , 

sin 2 z/ sinh 2 /i [A A P A (cosh /x) + P A Q A (cosh/x)] * (B.16) 
[C x P^cosv)+D x Q x (co S v)\ A^O . 

Pa and Q x are again solutions of the Legendre equation. Here, with reflection symmetry 
and integer A ()B.16|) looks like 

{Aq cosh /x + B A = , 

sin 2 v sinh 2 /xP^cos z/)L4 2m P2 m (cosh/x) (B.17) 
+P 2m Q' 2m (cosh /x)] A = 2/ + l . 

Oblate spheroidal coordinates : 

R = Ro cosh /x cos v , Q = R sinh /x sin v (B.18) 
cosh /x cos v jV M(U + cosh /x cos z/ — sinh /x cos z/iV^ + cosh /x sin z^iVj, = (B.19) 

Here, a separation with A(A — 1) G C, 3?A > —1 leads to 

{[A sinh /x + P ] [Co sin v + P ] A = , 

cos 2 v cosh 2 /x[A A P A (i sinh /x) + B x Q' x (i sinh/x)]* (B.20) 
[C A P^sinz/) + P A Q A (sinz/)] A^O . 

where, once again, P A and Q\ a Legendre functions of the first and second kind. By 
assuming integer A and reflection symmetry, ()B.20j) can be written as 

A sinh /x + Po A = , 

\l/ A (/x, v) = ^ cos 2 v cosh 2 /xP 2Z+1 (sin z/) L4 2 z + iP 2m (i sinh/x) (B-21) 

+P 2 m<2 2 m( isin M] A = 2/ + l. 
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